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POPULAR SUMMARY 

Every week, without realizing it, the U.S. population carries out a climate-change 
experiment by modifying its activities depending on the day of the week. It is well 
documented, for instance, that aerosol pollution generally decreases on weekends, 
apparently associated with reduced truck transport and power generation. Most urban 
dwellers are aware that pollution levels tend to be higher on weekdays, and monitoring 
stations maintained by the Environmental Protection Agency (EPA) at many locations 
around the U.S. show this. It is also now well established that aerosols affect the 
development of clouds. It is therefore natural to ask whether, on average, rainfall 
changes with the day of the week as a result of the experiment. Rain is highly erratic in 
its behavior, however, and it's to be expected that answering this question will require a 
considerable amount of observation to establish a pattern. 

The Tropical Rainfall Measuring Mission (TRMM) satellite has been orbiting the earth 
since the end of 1997. It is unique in carrying not only a meteorological radar capable of 
detecting rain from space, but other rain-sensitive instruments as well. Because the 
TRMM instruments have been observing the earth for over 10 years, and its 
measurements of rain are unusually trustworthy, the data it is supplying have become a 
gold mine for studies of changes in rainfall patterns. 

When we looked at eight years (1998-2005) of daily averages by TRMM of summertime 



rain over the southeast U.S., we found that more rain, and more intense rain, occurs in the 
middle of the week than on weekends in this region. Simultaneous reduction in rainfall 
over the nearby Atlantic and Gulf of Mexico is also seen. Changes are especially 
dramatic for afternoon rain (noon to midnight), and are statistically highly significant. 

We found that average wind patterns in these areas also change with the day of the week, 
consistent with the idea that the extra storm activity over the SE U.S. pulls more air in 
near the surface, carries it upwards, then pushes it outward and away at higher altitudes. 
We believe this outflow of air, squeezed dry by the storms, descends over the nearby 
waters and is responsible for the suppression of rain that we see occurring there at about 
the same time that rain is increasing over the SE U.S. The TRMM radar provides a 
different kind of evidence: it shows that storms tend to rise to higher altitudes in the 
middle of the week, and that rain occurs over larger areas. The changes in daily total 
rainfall also show up in rain-gauge measurements over SE U.S., though not as strongly as 
what is detected by the satellite instruments. 

After examining possible explanations, we have concluded that these changes are almost 
certainly due to higher midweek pollution levels. The observed rain intensification over 
land appears to confirm recent theories that pollution aerosols can cause summertime 
storms to grow more vigorously and possibly trigger development of nearby storms. 
Aerosols cause cloud drops to be smaller, which are able to rise higher where they freeze 
and release energy in the process, providing a kind of "after-burner" effect that energizes 
storm growth above what would occur without the aerosols. This effect is seen in 
numerical models of storm development, but it needs hot moist environments like the 
summertime SE U.S. 

Interestingly enough, rain gauge and wind analysis data suggest that the weekly cycle 
became strong enough to see only during the past few decades. Other research has 
suggested that heavy rainfall events have been occurring more often in recent decades. 
The TRMM observations of weekly variations in rainfall reported here suggest that some 
of this decadal increase may be due to increasing aerosol pollution levels during the past 
century rather than the changes in gas concentrations believed to be causing global 
warming. 
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ABSTRACT 

Tropical Rainfall Measuring Mission (TRMM) satellite data show a significant midweek 
increase in summertime rainfall over the southeast U.S., due to afternoon intensification. 
TRMM radar data show a significant midweek increase in rain area and in the heights 
reached by afternoon storms. Weekly variations in model-reanalysis wind patterns over 
the region and in rain-gauge data are consistent with the satellite data. A midweek 
decrease of rainfall over the nearby Atlantic is also seen. EPA measurements of 
particulate concentrations show a midweek peak over much of the U.S. These 
observations are consistent with the theory that anthropogenic air pollution suppresses 
cloud-drop coalescence and early rainout during the growth of thunderstorms over land, 
allowing more water to be carried above the 0°C isotherm, where freezing yields 
additional latent heat, invigorating the storms— most dramatically evidenced by the shift 
in the midweek distribution of afternoon-storm heights— and producing large ice 
hydrometeors. The enhanced convection induces regional convergence, uplifting and an 
overall increase of rainfall. Compensating downward air motion suppresses convection 
over the adjacent ocean areas. Pre-TRMM-era data suggest that the weekly cycle only 
became strong enough to be detectable beginning in the 1980's. Rain-gauge data also 
suggest that a weekly cycle may have been detectable in the 1940's, but with peak rainfall 
on Sunday or Monday, possibly explained by the difference in composition of aerosol 
pollution at that time. This "weekend effect" may thus offer climate researchers an 
opportunity to study the regional climate-scale impact of aerosols on storm development 
and monsoon-like circulation. 
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Abstract. Tropical Rainfall Measuring Mission (TRMM) satellite data show a signif- 
icant midweek increase in summertime rainfall over the southeast U.S., due to afternoon 
intensification. TRMM radar data show a significant midweek increase in rain area and 
in the heights reached by afternoon storms. Weekly variations in model-reanalysis wind 
patterns over the region and in rain-gauge data are consistent with the satellite data. 

A midweek decrease of rainfall over the nearby Atlantic is also seen. EPA measurements 
of particulate concentrations show a midweek peak over much of the U.S. These obser- 
vations are consistent with the theory that anthropogenic air pollution suppresses cloud- 
drop coalescence and early rainout during the growth of thunderstorms over land, al- 
lowing more water to be carried above the 0°C isotherm, where freezing yields additional 
latent heat, invigorating the storms — most dramatically evidenced by the shift in the mid- 
week distribution of afternoon-storm heights — and producing large ice hydrometeors. The 
enhanced convection induces regional convergence, uplifting and an overall increase of 
rainfall. Compensating downward air motion suppresses convection over the adjacent ocean 
areas. Pre-TRMM-era data suggest that the weekly cycle only became strong enough to 
be detectable beginning in the 1980’s. Rain-gauge data also suggest that a weekly cy- 
cle may have been detectable in the 1940’s, but with peak rainfall on Sunday or Mon- 
day, possibly explained by the difference in composition of aerosol pollution at that time. 

This “weekend effect” may thus offer climate researchers an opportunity to study the 
regional climate-scale impact of aerosols on storm development and monsoon-like circu- 
lation. 


1. Introduction 

The effect of pollution on rainfall has been observed 
to depend both on the type of pollution and the pre- 
cipitating environment [ Rosenfeld , 1999, 2000; Phillips 
et al 2002; Jacobson and Kaufman , 2006]. Pollution 
aerosols have been documented to suppress precipita- 
tion from shallow clouds (cloud heights below about the 
— 10° C isotherm) [ Albrecht , 1989; Rosenfeld , 1999, 2000; 
Rosenfeld et al. , 2002]. When polluted clouds develop to 
greater heights, however, as often happens in the sum- 
mertime over land, Rosenfeld suggested [Williams et al ., 
2002; Andreae et al ., 2004] that suppressed rainout en- 
ables unprecipitated water to reach greater heights where 
freezing can release additional latent heat and further in- 
vigorate the cloud updrafts. This might in turn delay 
the onset of precipitation and the development of down- 
drafts and so prolong the growth of the convective cloud, 
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allowing more water vapor to be ingested and further 
invigorate the storms [Rosenfeld, 2006]. Recent cloud 
simulations by Khain et al. [2005], Seifert and Beheng 
[2005], Wang [2005], Lynn et al. [2005], van den Heever 
et al. [2006], and Teller and Levin [2006] support this sug- 
gestion: in a moist, unstable atmosphere such as prevails 
during the summer in the southeast (SE) U.S., pollution 
aerosols can induce clouds to develop stronger updrafts 
and downdrafts, grow taller, trigger secondary storm de- 
velopment, and produce more rain. This is supported 
by satellite observations showing that clouds develop to 
greater heights in more polluted air masses [Koren et al . , 
2005]. 

It is well established [Simmonds and Keay , 1997; Sal- 
cedo et al. , 1999; Marr and Harley , 2002; Beirle et al ., 
2003; Blanchard and Tanenbaum , 2003; Bae et al ., 2004; 
Jin et al. , 2005] that pollution levels change with the day 
of the week in many urban areas, generally attributed to 
changes in vehicular traffic, though variations in power 
generation may also play a role. We provide examples 
of the widespread nature of this weekly variation over 
the U.S. in the next section. Such variations serve, in 
effect, as repeated experiments on the consequences of 
pollution. 

Weekly variations in temperature, pressure, cloud 
characteristics, hail and lightning are observed in many 
areas [ Lawrence , 1971; Fujibe, 1987; Gordon , 1994; Cer- 
veny and Balling , 1998; Dessens et al ., 2001; Marr and 
Harley , 2002; de F. Forster and Solomon , 2003; Mulla- 
yarov et al., 2005; Jin et al., 2005; Gong et al., 2006; 
Baumer and Vogel , 2007]. Searches for a weekly cycle in 
urban precipitation, however, have yielded mixed results 
[Cehak, 1982; Simmonds and Kaval, 1986; Simmonds and 
Keay, 1997; DeLisi and Cope, 2001; Jin et al, 2005; Gong 
et al. , 2006] . Over the Atlantic near the east coast of the 
U.S., Cerveny and Balling [1998] (“CB” hereafter), using 
rain estimates from the Microwave Sounding Unit (MSU) 
on TIROS-N satellites for 1979-1995, found a weekly cy- 
cle in precipitation with peak rainfall occurring on Sat- 
urday and minimum rainfall on Monday. A recent paper 
by Baumer and Vogel [2007] finds a weekly cycle in a 
number of meteorological variables measured at 12 sta- 
tions in Germany. Precipitation averaged over 15 years 
(all months) seems to peak on Saturday, as does cloud 
amount, but the daily averages used to find the peaks 
are somewhat noisy. 

Here we examine rainfall statistics from the Tropi- 
cal Rainfall Measuring Mission (TRMM) satellite [Kum- 
merow et al, 2000] over the southern U.S. and adjacent 
waters. The TRMM satellite has been orbiting the Earth 
since late 1997. It is unique in that it carries a meteoro- 
logical radar that can be used to improve the rain esti- 
mates made with its passive microwave instrument, the 
TRMM Microwave Imager (TMI). TRMM’s orbital plane 
is inclined 35° with respect to the Equator. The TMI can- 
not see poleward of 40° . Because it is in a low-inclination 
orbit, however, it views areas between 40 S and 40 N at 
all hours of the day over the course of 46 days, and is 
able to build up a statistical picture of changes in rain- 
fall statistics both with the time of day and season. 

We find a moderately significant midweek increase in 
daily TMI precipitation estimates over a large portion 
of the SE U.S. during the summer, and a highly signifi- 
cant midweek increase in afternoon rainfall over the area. 
TRMM radar measurements of storm heights (echo-top 
heights) over this area show a substantial midweek in- 
crease consistent with the intensification hypothesis out- 
^ lined above, and reanalysis winds show a midweek in- 
9 crease in low-level wind convergence, mid-level vertical 
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wind velocity, and upper-level wind divergence over the 
area, consistent with stronger midweek convection. Rain 
data from gauges and model reanalysis are also consistent 
with the satellite observations. Over the nearby waters, 
the reverse is found: a highly significant weekend peak in 
rainfall over the Atlantic and a weaker one over the Gulf 
of Mexico. 

These observations open a window to the effects of an- 
thropogenic pollution on regional rainfall at climatologi- 
cal scales. The intensification over land and what appear 
to be compensating effects over adjacent waters support 
the suggestion that air pollution suppresses cloud drop 
coalescence and early rainout, allowing subsequent invig- 
oration of the thunderstorms and outflow aloft. It sug- 
gests that summer rainfall on large scales both increases 
and intensifies as pollution levels rise. The increases in 
frequency of heavy rainfall events over the U.S. during 
recent decades found by Groisman et al. [2004] may be 
partly explained by this mechanism. 

In the following section we analyze data for particulate 
aerosols and show that the aerosol concentration gener- 
ally peaks in the middle of the week over large areas of 
the U.S. In Section 3 TRMM data are used to look for 
a weekly cycle in precipitation, and a distinct weekly cy- 
cle is found over the SE U.S. and neighboring waters. 

When rainfall is separated into morning and afternoon 
amounts, the weekly changes are found to be statistically 
highly significant. In Section 4 we examine TRMM radar 
data for echo-top heights and find a dramatic midweek 
increase in rain frequency and heights reached by the 
most intense storms. In Section 5 we find that reanalysis 
data for atmospheric winds indicate weekly variation in 
wind patterns over the same areas that is consistent with 
the weekly cycle in convective activity. We also exam- 
ine rain-gauge data and reanalysis estimates of rainfall * 
and compare them with TRMM. In Section 6 we look for 
signs of a weekly cycle in daily rain estimates from gauges 
and reanalysis data prior to the TRMM era, Results are 
discussed in Section 7, and our conclusions presented in 
the final section. Two appendices provide additional in- 
formation concerning the analysis. 

2. Weekly Cycle of PM i0 and PM 2 . 5 

Pollution varies with the day of the week in many areas 
of the world. Some examples of observations of a weekly 
cycle for various types of pollution are provided, for in- 
stance, by Cleveland et al. [1974], Simmonds and Keay 
[1997], Cerveny and Balling [1998], Salcedo et al [1999], 

Diem [2000], Marr and Harley [2002], Bae et al. 2004 , 

Harley et al. [2005], Jin et al. [2005], and Shutters and 
Balling [2006] . Beirle et al. [2003] and Beirle et al. [2004] 
provide a nice global perspective using satellite observa- 
tions of the weekly cycle in NO 2 production. 

The widespread nature of the weekly cycle in pollution 
over the U.S. is evident in an analysis of data for par- 
ticulate concentrations obtained from the Environmen- 
tal Protection Agency (EPA) Technology Transfer Net- 
work Air Quality System [USEPA, 2006]. EPA measure- 
ments of concentrations of particulates smaller than 2.5 
jam (PM2.5) and 10 jam (PM10) have been made at many 
locations around the U.S. at regular time intervals rang- 
ing from days to hours. They are generally reported in 
units of fag m -3 . Discussions of PM2.5 and PM10 mea- 
surements may be found in Wilson [2002] and Dye et al. 

[2003]. 

EPA data for years 1998-2005, June- August, were an- 
alyzed. They were current as of 24 Aug 2006. Where 
available, daily values of PM2.5 and PM10 were used. In 
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addition, for those sites where hourly rather than daily 
values of PM 10 were available, the hourly values were av- 
eraged to daily values, excluding days missing more than 
6 hourly values, as well as any days containing negative 
hourly readings during that day. Site data for years with 
fewer than 69 days (75% of 92 days) were discarded. With 
these criteria, 212 sites provided 1-8 summers of data for 
PM 2 . 5 , while 386 sites provided data for PM 10 . Some 
“sites” are quite close to each other physically and are 
mostly distinguished by the kind of measuring appara- 
tus used. The number of summers of data available at 
collocated sites may also differ. 

Daily data for each of the sites were fit to 7-day si- 
nusoids. Estimates of the statistical significance level p 
of the sinusoidal amplitudes were obtained following the 
technique described by Bell and Reid [1993]. Some de- 
tails of the techniques are given in the next section and in 
Appendix A. The interpretation of a significance level p 
for a given site is that a spurious weekly-cycle amplitude 
as large as what is observed, under the hypothesis that 
mo real cycle is present, could occur “by accident” with a 
probability p determined from an estimate of the level of 
random Variability in weekly excursions. Smaller values 
of p suggest that the cycle is more likely to be “real” . 

The maps in Figures 1 and 2 show the results of an- 
alyzing the EPA data. Each vector originates from the 
location of a monitoring site, with its direction indicating 
the day of the week on which the 7-day sinusoidal fit to 
the daily averages peaks (see key to directions on right 
side of maps) , and with the lengths and colors of the vec- 
tors indicating the statistical significance of the weekly 
cycle. Vector lengths are proportional to (— lnp) 1/2 , but 
are not allowed to exceed the length corresponding to 
p = 0.01. A key to the length scale is provided at the 
bottom right corner of the maps. 

It is clear from Figure 1 that there is a widespread 
tendency for the average concentration of PM 10 to peak 
Tue-Thu, although there are almost certainly real differ- 
ences in the strengths and phases of the cycles depending 
on location. A midweek peak in PM2.5 is also generally 
evident over much of the U.S. (Figure 2), with notable 
anomalies in the center of the country. The weekly cycle 
in PM2.5 appears to be weaker than in PM10, perhaps 
because lighter particles settle out more slowly and are 
more affected by weather, so that their concentration is 
less strongly tied to the weekly cycle in emissions. (The 
lengths of the time series at PM2.5 and PM10 sites are 
similar, so the higher level of “noise” in the weekly cy- 
cle of PM2.5 is not due to differences in the amount of 
data.) Based on the amplitudes of fits of the data to a 
sinusoid peaking on Wednesdays, we find that the ampli- 
tudes of the weekly cycles of PM 10 are typically of order 
±10% of the mean concentrations, while those of PM2.5 
are of order ±5% of the mean; but, again, both means 
and amplitudes surely vary with location. 

A number of studies have tried to characterize the 
weekly changes in particulate concentrations. Bae et al. 
[2004], for example, find that elemental carbon particu- 
late concentrations are about 20% lower on weekends in 
the neighborhood of East St. Louis, IL. Marr and Harley 
[2002] find that, in the Sacramento/ San Joaquin Valleys 
in California, heavy truck traffic, which mostly uses diesel 
fuel, is substantially smaller on weekends than weekdays, 
while daily passenger and fight- truck traffic do not vary so 
much. The results of Wahlin et al. [2001] suggest some- 
thing similar in the vicinity of Copenhagen, Denmark. 
Charron and Harrison [2005] find that heavy-duty diesel 
vehicles are responsible for substantial weekly variations 
in PM2.5 and PM10 on a major arterial road near Lon- 
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don, Great Britain. Perry and Owens [2001] document 
weekend minima in traffic levels, CO and nitrogen oxide 
concentrations, and power generation in the Charlotte, 

North Carolina region. 

It is worth noting that the weekly cycle in total column 
aerosol amount may not be nearly as strong as the weekly 
cycle in near-surface aerosol concentrations as reflected in 
the measurements of PM2.5 and PM 10 described above. 

We examined satellite estimates of aerosol optical thick- 
ness (AOT) over an area of the southeast U.S. bounded 
by 32.5N-40.0N, 100W-80W from the Moderate Resolu- 
tion Imaging Spectroradiometer (MODIS) aboard both 
NASA’s Terra and Aqua satellites [Remer et al . , 2005], 
using “Collection 4” gridded products MOD08-D3 and 
MYD08JD3 downloaded using the GES-DISC Interactive 
Online Visualization ANd aNalysis Infrastructure (Gio- 
vanni) as part of the NASA’s Goddard Earth Sciences 
(GES) Data and Information Services Center (DISC). 

We found that the weekly cycles in both satellites’ day- 
time estimates of AOT were not statistically significant, 
though sinusoidal fits to both satellite estimates peaked 
on Wednesday. (The nominal overflight time of Terra is 
10:30 AM and that of Aqua is 1:30 PM.) 

3. Weekly Cycle in TRMM Rainfall 

To look for signs of the influence of weekly variations 
in human activity on precipitation, averages for each 
day of the week based on 8 years (1998-2005) of data 
for summertime (June-August) rain rates, estimated us- 
ing version 6 of the TMI retrieval algorithm [Kummerow 
et al., 2001; Olson et a/., 2006], were obtained for each 
2.5° x 2.5° grid box viewed by TRMM in the vicinity of 
the continental U.S. (Some results for 2006 will be dis- 
cussed later.) 

The daily averages r(t) for each grid box were fit to a 
sinusoidal form 

r(t) = r 0 + r 7 cos [u> 7 (t - </> 7)] (1) 

with uo 7 = 27r/(7days), t the time measured in days, ro 
the average rain rate, r 7 the amplitude of the weekly cy- 
cle, and (j ) 7 the day of the week when the sinusoidal fit 
peaks. The statistical significance of the amplitude r 7 
under the hypothesis that there is no weekly cycle (r 7 = 

0) was obtained using a technique described in Appendix 
A. Figure 3a shows a map of the phase cj ) 7 obtained using 
this fitting procedure. The strongest colors (topmost in 
color bar) indicate significance levels of r 7 with p = 0.05 
or better (which it may be helpful to think of as equiv- 
alent to amplitudes larger than “2 sigma” for normally 
distributed errors). Although there is considerable vari- 
ability in the phase (much of which might be explained 
by sampling error due to the length of the dataset), it 
appears that there is a tendency for average rain rates to 
peak during the middle of the week (Tue-Thu, predom- 
inantly reddish hues) over the continental U.S., and to 
peak Sat-Mon over the nearby Atlantic (predominantly 
bluish hues) and perhaps the Gulf of Mexico. There also 
appears to be a tendency for the weekly cycle to weaken 
near the coasts, in the sense that significance levels tend 
to be low there, or for the peaks to shift to Friday. Re- 
sults of studies of data from coastal cities mentioned ear- 
lier looking for a weekly cycle in precipitation may have 
been inconclusive for this reason. 

Even with 8 years of data, the statistical uncertainty 
in the weekly cycle in many of the grid boxes is large. 

We therefore try to increase the signal-to-noise ratio by 
averaging over larger areas, guided by our physical under- 
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standing of where and how the weekly cycle in pollution 
is likely to affect precipitation. Summer precipitation 
over the U.S. differs markedly in the eastern and western 
halves, due in part to the effects of the moisture brought 
to the eastern half from the Gulf of Mexico by prevailing 
winds [see Jin et a/., 2005, for example] and in part to 
the differences in topography. We have therefore calcu- 
lated average rain rates for the five averaging areas A-E 
shown in Figure 3b. Grid boxes containing substantial 
amounts of coastline are excluded, partly based on an 
apparent reversal of the phase at the coasts and partly 
based on concerns about the change in the TMI retrieval 
algorithm that occurs at the coast [e.g., Kummerow et al . , 
2001]. Area C is identical to the region examined by CB 
except for two of their 2.5° grid boxes too far north to 
be seen by TRMM. 

Figure 4a shows the average rain rates for each day of 
the week for areas A-C. Averages are obtained by sum- 
ming rain- rate estimates for all TMI footprints falling 
within the area during the averaging period, and divid- 
ing by the number of footprints contributing to the sum. 
(One-sigma error bars for Area-B averages are estimated 
using the resampling technique described in Appendix 
A.) The SE-U.S. (area-B) average daily rain rates for 
Tue-Thu are higher than for Sat-Mon, with a maximum 
on Tue. The SW U.S. (area A) may have a tendency to 
have higher midweek rain rates, but the signal is very 
weak. In contrast, the coastal Atlantic (area C) shows 
strongest rain rates Sat-Mon and lower rain rates Tue- 
Thu — almost exactly opposite to what is happening over 
land. The behavior over area C is similar to what CB 
found, but their results were based on a different sensor 
for an earlier time period (1979-1995), and they averaged 
data over all seasons whereas our averages include only 
summers. [A fit of Eq. (1) to 7 complete years of TMI 
data (1/98-12/04) for area C peaks on Sunday but has 
little statistical significance (p = 0.75). Our dataset does 
not show a statistically significant all-season weekly cy- 
cle analogous to what CB found, possibly because it is 
shorter in length than theirs.] 

Figure 5a shows the average rain rates for areas C- 
E. Statistics for area C have been repeated there to aid 
comparisons. Both the Atlantic region D (east of region 
C) and the Gulf of Mexico (area E) seem to show weekend 
maxima similar to that of C, but weaker. 

3.1. Weekly cycle in rainy area and intensity 

The weekly changes in rain rate are due in part to 
changes in the average fraction of the area covered by 
rain and in part to changes in the intensity of rain where 
it is raining. This can be seen in Figures 4b, c and 5b, c. 
The average fraction of area with rain is represented by 
the daily average fraction of TMI footprints (nominally 
of order 10 km in diameter) with detectable rain in them 
for each area, while the average intensity of rain is rep- 
resented by the ratio of the mean rain rate (Figures 4a 
and 5a) to the mean fraction- with- rain (Figures 4b and 
5b) , equivalent to the rain rate averaged over rainy areas 
only. Linear regression of changes in average rain rate 
for each day against changes in rainy area and changes 
in intensity suggests that increases in average rain rate 
are about 1/3 due to increases in areal coverage by rain 
and 2/3 to rain intensity. Data for TRMM Precipita- 
tion Radar (PR) storm heights that will be discussed in 
Section 4 also indicate both a midweek increase in storm 
intensity, as inferred from increased storm heights, and 
increased areal coverage by rain. 


Figure 3b 


Figure 4a 


Figure 5a 


Figures 4b, c and 5b, c 
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3.2. Statistical significance of cycles 

The question immediately arises whether these cy- 
cles are artifacts of the observational characteristics of 
TRMM, or accidents arising from the natural variability 
of precipitation. Evidence that the sampling pattern of 
TRMM is unlikely to have generated spurious weekly cy- 
cles of the magnitude seen here is described in Appendix 
B. The possibility that the weekly oscillations seen in 
Figs. 4 and 5 are simply residues of random rain events 
with no favored days of the week is evaluated using a 
resampling (bootstrap) technique that attempts to take 
into account the multi- day time correlations of rain and 
the special sampling pattern of the TRMM satellite. It 
is described in Appendix A. Figure 6 displays the re- 
sampling results for the 5 areas as a “clock plot.” The 
phases 07 of the fits are indicated by angular position, 
whose meaning is given by the labeling along the clock 
perimeter. The statistical significance p of each cycle, es- 
timated using resampling, is indicated by distance from 
origin, and is proportional to (— 1 np) 1//2 . Circles with 
radii corresponding to various values of p are drawn to 
define the radial scale. 

The weekly cycle of the average over the SE U.S. (area 
B) is significant at the p = 0.10 level. (Parameter values 
of fits for each area are given in Table 1.) The coastal At- 
lantic (area C) cycle is significant at the p = 0.035 level. 
The other two areas over water, D and E, have weekly 
amplitudes that are above “1-sigma” in size. A weekly 
cycle over the SW U.S. (area A), if present, is not strong 
enough to be detected with any confidence. In fact, a 
weakened response over the SW U.S. (area A) is consis- 
tent with simulations by Khain et al [2005], which show 
that invigoration of convection by pollution is diminished 
in drier conditions, and may even reverse. 

[It may be helpful to the reader to be reminded that 
plots like Figure 6 can serve two purposes: The length 
of a vector gives the significance level p of the amplitude 
f7, indicated by the circles labeled by various values of 
p, and serves to test the hypothesis 7*7 = 0. Once one 
has accepted the hypothesis that there is a weekly cycle 
(7*7 7^ 0) and one wants confidence limits for 7*7, 07, the 
same plot can provide this if one draws circles centered on 
the end of the vector with radii equal to those of the cir- 
cles labeled by p. Thus, the three circles in Figure 6 can 
be used to provide the 70%, 90%, and 97% confidence- 
limit circles for the 5 vectors. It should be remembered 
that each vector rj is plotted scaled by its own “noise” 
estimate <77 (see Appendix A), and so quantitative com- 
parison of two amplitudes 7*7 requires replotting of the 
data with a uniform scale. Additional discussion may be 
found, for example, in Collier and Bowman [2004].] 

3.3. Morning vs. afternoon statistics 

The significance level of the SE-U.S. weekly cycle does 
not quite reach the “canonical” p = 0.05 level. We note, 
however, that the explanation for the intensification of 
storms by pollution proposed here would suggest that 
the intensification should be greatest when atmospheric 
instability is greatest, since moist parcels of air are car- 
ried highest in such an environment. This suggests that 
the intensification should be greatest for afternoon rain- 
fall, and this is indeed what is seen. 

Figure 7 shows the average rain rate for each day of the 
week for morning rain (0000-1200 LT) and afternoon rain 
(1200-2400 LT) over area B. Local time (LT) is computed 
from the longitude A (—180° < A < 180°) of the centroid 
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of each area as 


LT = UT + A/(15°/h) , (2) 

where UT is Universal Time. The weekly cycle of af- 
ternoon rain is quite strong. It peaks on Tuesdays, and 
has a very high significance level p = 0.0012 when esti- 
mated using resampling. The afternoon weekly ampli- 
tude is about twice as large as the morning amplitude, 
which has a significance level p = 0.045 and peaks on Sat- 
urdays. The impact of the afternoon cycle on daily totals 
thus seems to be somewhat mitigated by compensating 
changes in morning rainfall, so that the weekly cycle for 
daily total rainfall is weaker than for the separate morn- 
ing/afternoon components. Ways in which precipitation 
during one part of the day can influence precipitation at 
other times, thereby affecting the overall diurnal cycle, 
are discussed by Betts and Jakob [2002]. 

Although not shown here, the relative contributions to 
the weekly change in afternoon rainfall from changes in 
rainy area and intensity are rather different from what 
was inferred from the 24-h totals in Section 3.1: we find 
about 75% of the weekly peak in rain rates is due to in- 
crease in rainy area and 25% to increase in rain intensity. 
Note, however, that over land the TMI estimates of peak 
intensities may be underestimated due to saturation of 
the 85- GHz channel due to the large amount of ice aloft. 

. The afternoon cycle is so strong that, when individual 
summers of data are fit to sinusoids, 5 of the 8 summers 
(1998-2005) peak on Wed, and the remaining 3 on Mon- 
Tue, as shown in Figure 8a. The sinusoidal fit to TMI 
data for year 2006, however, peaks on Saturday, with a 
significance level p = 0.08 estimated using resampling. 
In examining the weekly cycles of afternoon TMI data 
for single summers, however, we should remember that 
because the TRMM satellite views the U.S. during after- 
noon hours for about 3 weeks and then spends another 
3 weeks observing during morning hours, only about 6- 
7 weeks of afternoon data are collected each summer. 
The estimates for individual years shown in Figure 8a 
are therefore not based on many weeks of data. To see 
statistical evidence of a weekly cycle, nature must pro- 
vide us with a set of trials that are evenly distributed over 
the days of the week, with similar conditions for storm 
development but varying pollution levels. If favorable 
conditions for rain occur at best a few times per week, 
a single summer of afternoon TMI observations is quite 
likely to fail to provide us with a set of “experiments” 
distributed evenly enough for a convincing weekly cycle 
to be visible in the data. 

To get a better feeling for the stability of the yearly 
estimates in Figure 8a, data from the TRMM Multi- 
satellite Precipitation Analysis (TMPA) product were 
analyzed. This TRMM product, 3B42 in the standard 
TRMM product list, incorporates both data from var- 
ious satellites carrying microwave instruments, includ- 
ing TRMM, and rain estimates using geosynchronous 
infrared (IR) data, tuned to match the microwave es- 
timates. The estimates are adjusted where possible to 
agree with monthly rain-gauge totals. The product pro- 
vides gridded rain fields every 3 hours, and is described 
by Huffman et al [2007]. The data were acquired us- 
ing the Giovanni interface mentioned in Section 2. The 
TMPA product increases the effective number of samples 
per summer by perhaps a factor of 4 or more over the 
TMI data alone. Over land, the accuracy of individual 
rain estimates is likely to be roughly comparable to that 
of TMI estimates. The accuracy of the IR estimates is 
the most problematic, but, on a global basis, the fraction 
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of the dataset for which IR estimates are the only infor- 
mation declines from about 55% in 1998 to 20% in 2004 
[G. Huffman , 2007 , private communication]. Although 
this dataset brings along with it a set of remote-sensing 
and algorithmic issues that are not easy to quantify, it 
provides a quasi-independent satellite perspective that 
we believe may be helpful. We have used the 3-hourly 
TMPA rain estimates averaged over a rectangular area 
that is close to that of area B, covering 32.5N-40.0N, 
100W-80W, obtaining fits to each summer of data for 
the available years 1998-2006. The results are plotted in 
Figure 8b. The estimates for the two most anomalous 
years 2005 and 2006 have moved closer to those for the 
earlier years and are weaker in strength. Such changes 
in the strength of the weekly cycle from year to year are 
consistent with variations we shall discuss later for data 
from rain gauges, and suggest that the weekly-cycle signal 
waxes and wanes in strength from summer to summer. 

(This is the only section in which data for 2006 is dis- 
cussed and the TMPA product is used.) 

3.4. Weekly and diurnal cycles in TMI averages 

A more detailed picture of the day-night differences 
in the weekly cycle can be seen in Figure 9, which shows 
how the weekly cycle of rain in areas B and C varies 
with the hour of the day. Note how the cycle over land 
(area B) is mirrored over the nearby ocean (area C). The 
afternoon maximum over land on Tuesday is paired with 
a minimum over the ocean with a delay of ~2 hours. Note 
also how much weaker the diurnal cycle is on Saturdays 
than on Tuesdays over area B. 

4. TRMM Radar Storm Heights 

The mechanism for the midweek intensification of 
storms proposed here suggests that more storms should 
reach higher altitudes Tue-Thu than Sat-Mon. The 
TRMM radar product 2A23 [TRMM PR Team , 2005] 
reports “storm height” for each radar observation when 
the PR algorithm determines that there is rain detected 
within the radar beam with a high degree of confidence. 
“Storm height” is the height of the highest point in the 
radar beam with detectable returns, measured relative to 
mean sea level. The PR beam is roughly 4-5 km in di- 
ameter. A detectable radar return indicates the presence 
of large water droplets or ice aggregates that have been 
carried aloft by strong vertical winds; they will begin to 
fall when the winds are no longer present. These storm 
heights are distinct from cloud- top heights, which might 
be much higher. There are a number of subtle issues 
involved in the interpretation of PR storm heights, due 
to effects such as measurement noise and horizontal tilt- 
ing of storms, for example, but since we are comparing 
statistics from different days of the week, these issues are 
unlikely to change the general conclusions we draw from 
these statistics. 

We collected PR storm-height statistics over the SE 
U.S. (area B) for JJA 1998-2005. “Weekend” (Sat-Mon) 
and “midweek” (Tue-Thu) statistics were obtained sep- 
arately for morning and afternoon time periods. The 
number of weekend and midweek PR footprints with 
valid data (including those with no detectable rain) was 
6.8 x 10 6 each for the morning hours, and 6.6 x 10 6 each 
for the afternoon hours. Table 2 gives the fraction of 
valid PR observations with rain in them for the various 
cases (i.e. , the average areal fraction with rain). It shows 
that the change in average rainy area from morning to af- 
ternoon (an aspect of the diurnal cycle of precipitation) 
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is much smaller on weekends than midweek, consistent 
with the TMI statistics in Figure 9, which shows that 
the diurnal cycle is weaker on weekends. It also shows 
that in the afternoon hours rain is detected about 40% 
more midweek than on weekends, while in the morning 
hours rain is detected 20% less. 

Figure 10a shows the frequency distributions of storm 
heights seen by the PR (number of storm heights in 1- 
km bins expressed as a fraction of the total number of 
PR footprints with detectable rain), for “weekends” and 
“midweek” during morning hours. The difference of the 
weekend distribution from the midweek distribution is in- 
dicated by the dashed curve. It shows that the distribu- 
tion of midweek storm heights shifts to higher altitudes 
relative to weekend heights. Figure 10b shows that a 
similar but much larger change in the distribution occurs 
during the afternoon hours. 

We can calculate the frequency with which storm 
heights exceed a given altitude a (i.e., cumulative dis- 
tribution) by summing the corresponding distribution in 
Figures 10a, b from the highest altitude down to a. The 
ratio of the midweek cumulative distribution to the week- 
end distribution then gives the relative probability that 
midweek storm heights will exceed a given altitude com- 
pared to weekend storm heights. The ratios are plotted 
for morning and afternoon storms in Figure 10c. The 
probability of afternoon storm heights exceeding 9 km is 
about 40% higher midweek than on weekends. In other 
words, not only is the average area covered by afternoon 
storms 40% higher Tue-Thu compared with Sat-Mon, 
but, once a storm begins to develop, it is 40% more 
likely to reach altitudes above 9 km on Tue-Thu than on 
Sat-Mon. We conjecture that the reason morning storms 
reach altitudes above 15 km relatively more often than 
afternoon storms is that the highest storms are also the 
biggest and longest-living, and the time for organization 
at these scales often means that they reach their peak 
altitudes after midnight (and are counted as “morning” 
storms in Figure 10). 

Deriving confidence limits for the frequency distribu- 
tions shown in Figures 10a, b is complicated by the spatial 
and temporal correlations of storm heights: the counts in 
a bin cannot be assumed to be independent, especially at 
lower altitudes. At the higher altitudes, however, where 
counts are infrequent, an assumption of independence of 
the individual counts is probably not so bad. If this as- 
sumption is made, error analysis of the distributions at 
altitudes above 9 km (not shown here) suggests that the 
change in the distributions from weekend to midweek is 
highly significant. The method for estimating the error 
bars in the ratios plotted- in Figure 10c is described in 
Appendix A, and should probably be considered trust- 
worthy only above 9 km. 

5. Cycles in Reanalysis and Gauge Data 

5.1. Reanalysis winds 

The weekly cycles in areas B and C are strong enough 
that it would be surprising if they were not reflected in 
corresponding variations in wind fields in the region. Us- 
ing version R-2 of the National Centers for Environmental 
Prediction/Department of Energy (NCEP/DOE) reanal- 
ysis data [Kanamitsu et al . , 2002] (referred to here as 
NCEP2), we obtained wind statistics over areas B and C 
that seem likely to reflect changes in convective activity 
in the areas: average surface (1000-hPa) convergence of 
the horizontal winds, vertical wind at 500 hPa, and wind 
divergence at 300 hPa. Note that the daily reanalysis 
data we use are averages from 00Z to 24Z. At typical 
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eastern U.S. longitudes, the data therefore represent av- 
erages over roughly 1800 LT the previous day to 1800 LT 
of the nominal (UT) day of the reanalysis. 

Figure 11 shows these reanalysis- wind statistics as a 
function of the day of the week, derived from summer 
(JJA) data for 1998-2005. The signal strengths of the 
weekly cycles over land are very high, as measured by 
the significance levels p (see figure caption); the cycles 
over the ocean are noisier. The weekly cycle in air mo- 
tion represented by these averages corresponds quite well 
with the behavior one would expect to accompany the 
vertical convection associated with the rainfall cycles ob- 
served by TRMM. Note that for area B the TMI rain rate 
and all but the 300-hPa divergence have a minimum on 
Sunday, while for area C all have a minimum on Tuesday, 
even though the satellite data and the NCEP-reanalysis 
data use very different measurement systems with very 
different sampling characteristics. 

5.2. Reanalysis and gauge rainfall 

Rain gauges measure rainfall more directly than any 
remote sensing technique, but only over an area the 
size of a dinner plate. Since variations in rain rate are 
spatially correlated, however, measurements by a single 
gauge reflect changes in rainfall over an area surround- 
ing the gauge, and so a sufficiently dense array of gauges 
can detect changes in spatially averaged rain rates. The 
National Climatic Data Center (NCDC) has developed 
a quality-controlled dataset of thousands of rain-gauge 
records across the U.S. [ Vose et al, 1992]. We used their 
product designated “Data Set 9300” from the Global His- 
torical Climatology Network (GHCN), version 1.0, up- 
dated through 2005. The dataset we used covers 1901- 
2005, though the number of gauges begins to diminish 
steadily as one goes back in time from 1950. [Version 2 
of this dataset was released after this analysis was com- 
pleted. See GHCN, 2007]. 

The NCEP2 reanalysis attempts to assimilate rain 
data by forcing soil moisture in the model to adjust to 
measured rain rates, supplied to the model as 5-day run- 
ning averages. The first version of the NCEP-reanalysis 
data [Kalnay et al , 1996], referred to here as NCEP1, dif- 
fers from NCEP2 in not attempting to assimilate rainfall 
data, but covers the longer period 1948-2005: NCEP2 
only covers 1979-2005. 

NCEP-reanalysis rain estimates were averaged over ar- 
eas B and C, and GHCN gauge data over area B, for 
1998-2005 summers. Figure 12a shows the mean rain 
rate for each day of the week for area B estimated by 
the TMI, NCEP reanalyses, and GHCN rain gauges (ap- 
proximately 1800 in number and fairly evenly distributed 
geographically). The different techniques clearly give dif- 
ferent mean rain rates, and the differences are probably 
due to real biases in the methods rather than simply due 
to sampling differences. All have a minimum on Sunday, 
but the peaks range from Tuesday to Thursday. Given 
the sampling errors in the various estimates represented 
by the error bars, it is possible that the differences in the 
TMI, reanalysis and gauge peaks may simply be due to 
sampling. Note that the weekly phase of the two ver- 
sions of the NCEP reanalysis seem to differ, even though 
they share a considerable amount of data, possibly due 
to the method of assimilating rain-rate data in NCEP2. 
Figure 12b shows the corresponding picture for area-C 
averages. Mean rain rates again differ among the differ- 
ent methods. Both NCEP-reanalysis daily means peak 
on Saturday, while the TMI peaks on Sunday. Table 3 
gives the parameters for fits to the various data and the 
significance levels p of the fits. 
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6. Historical Behavior of Weekly Cycle 

We have so far looked for weekly cycles in precipitation 
and other fields during the TRMM era, which started in 
December 1997. The NCEP1 reanalysis data extend back 
to 1948, however, and the GHCN rain-gauge data extends 
back to 1901, though coverage dwindles from about 1200 
gauges in the later 1940’s to about 400 in 1901. We can 
therefore look for a weekly cycle in these data prior to 
TRMM. 

Examining the historical behavior of the strength of 
the weekly cycle is subject to a number of uncertainties. 
From a physical perspective, pollution levels and types of 
pollution have changed with the decades due to changes 
in regulation, technology, and population. If it is indeed 
pollution causing most of the weekly cycle we see in the 
data, we cannot be sure the pollution types of years ago 
had the same impact as current pollution types are hav- 
ing. From a statistical perspective, the strength of the 
signal we have seen for daily averages in the 8- year period 
1998-2005 over area B varies with the quantity examined, 
with significance levels from p = 0.01 to p = 0.1 (see cap- 
tion to Figure 11 and Table 3). If such significance levels 
are characteristic of the signal in earlier years, then, for 
example, an examination of ten 8-year periods is almost 
guaranteed to produce cases with estimated significance 
levels of p = 0.1 even if they are entirely spurious. 

The TRMM data suggest that our best hope of learn- 
ing whether a weekly cycle has been present in ear- 
lier years would be to restrict our analysis to afternoon 
data, since the weekly cycle seems so much stronger in 
the afternoon data. The NCEP-reanalysis data are in 
fact available at 6-hourly intervals, seeming to offer us 
a chance to look at atmospheric behavior for afternoon 
hours alone. When we looked at NCEP-reanalysis data 
for 00Z (corresponding to “the afternoon,” about 1800 LT 
over the SE U.S.) and 12Z (corresponding to about 0600 
LT) results seemed generally to be consistent with what 
the TRMM data lead us to expect for the last decade, 
but are not so easy to interpret for earlier decades. Since 
the NCEP data assimilation relies heavily on radiosonde 
data that are supplied at best only twice per day, and the 
models themselves have some difficulty in representing 
the diurnal cycle of convective behavior [e.g., Janowiak 
et al . , 2007], it is possible that current data assimilation 
products are not yet suited to our needs. Rain-gauge 
data at hourly intervals may be able to help here, but re- 
quires analysis of an order of magnitude more data than 
we have looked at so far, and we have not done this yet. 

With these caveats in mind, we present a few examples 
of what we see for the historical behavior of the weekly 
cycle in rain-gauge and NCEP-reanalysis data. We use 
the signal-to-noise ratio rj/cri to represent the strength 
of the weekly cycle, with the quantity cr 7 defined in Ap- 
pendix A representing the estimated sampling variability 
in the weekly- amplitude parameter 77 , from which the 
significance level p can be obtained using the formula 
p — exp [— (77/(77) 2 ]. [A plot of <77 for the area-B rain- 
gauge average for each year (not shown) is fairly con- 
stant, and does not rise in the earlier years as might be 
expected if the decline in the number of gauges for years 
prior to 1950 caused the representativeness of the area 
average to diminish with the number of gauges.] 

Figure 13a shows the signal-to-noise ratio of weekly- 
cycle fits for daily area- averaged gauge data over area B, 
plotted as far back as 1908. The statistics of fits plotted 
at year y are based on the 8 years y — 7 to 7/, using only the 
summer data. We have chosen an 8-year window solely 
because it is the size of the window used in the TRMM 
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studies in the earlier sections of this paper. Different 
windows can produce significantly different-looking plots. 

(One example is shown later.) Two horizontal lines are 
drawn on the plot indicating the values of 7*7 jcti for which 
p — 0.5 and 0.05. If there were no weekly cycle, we would 
expect half the ratios to fall below p — 0.5 and half to 
fall above it. We would expect 5% of the ratios to exceed 
the p = 0.05 level by accident. The figure suggests that 
a weekly cycle might have begun increasing to detectable 
levels in the 1980’s, but this is far from clear. The day 
of the week <fi 7 when the weekly fit peaks is plotted in 
Figure 13b. Note that the phase of the fits becomes more 
erratic for low signal- to-noise ratios (e.g., year 2001 in 
Figure 13). 

There is a hint in Figure 13 that a weekly cycle might 
have appeared in the decades around 1940, peaking Sun- 
day/Monday, then diminished and reappeared in the late 
1970’s with a peak closer to the middle of the week. It 
is possible that this might be related to the changes in 
the composition of air pollution over these decades from 
absorbing to less absorbing aerosols: Lefohn et al. [1999] 
suggest that sulfur emissions in the U.S. might have sta- 
bilized beginning around 1920, possibly due to the impact 
of government regulatory action, while Novakov et al. 

[2003] suggest that black carbon emissions from diesel 
engines began to rise in the 1950’s while black carbon 
emissions from other sources decreased. Diesel particles 
are generally smaller than 100 nm. Their optical cross 
section is very small, but they are still very active as 
cloud condensation nuclei (CCN) when chemically ma- 
ture [Dusek et al . , 2006]. This suggests that there might 
have been a transition from large absorbing aerosols dur- 
ing the first half of the 20th century to less absorbing, 
smaller and more numerous particles towards the end of ' 
the century. Absorption by aerosols reduces solar heat- 
ing of the surface and hence suppresses convection [Koren 
et al ., 2004], while the CCN activity of aerosols enhances 
deep convection [ Rosenfeld , 2006]. This can potentially 
explain a transition from a weekly cycle that is dominated 
by aerosol absorption suppressing midweek convection in 
the middle of the 20th century to the period towards 
the end of the 20th century in which the elimination of 
visible black smoke causes a changeover to midweek in- 
tensification due to the increasing number of small, sub- 
visible CCN along with decreasing solar absorption by 
the aerosols [ Wild et al . , 2005] . 

Figure 14 shows the signal-to-noise ratios and phases 
for NCEP1 surface rainfall, surface wind convergence, 
300-hPa wind divergence, and 500-hPa vertical velocity 
with an 8-year running window. There appears to be 
a tendency for the weekly signal to begin strengthening 
in the 1980’s in all plots. This tendency appears much 
more dramatic when a window longer than the 8-year 
window used in Figures 13-14 is used: Figure 15 shows 
the signal-to-noise ratios for running 15-year fits to the 
gauge and NCEP1 reanalysis fields. It should be noted 
that the NCEP1 reanalysis incorporated TIROS (Televi- 
sion Infrared Observation Satellite) Operational Vertical 
Sounder (TOVS) data beginning in about 1979 [Kalnay 
et al . , 1996], and it is conceivable that this might have 
increased the signal-to-noise ratio of a weekly cycle rel- 
ative to what could be observed prior to 1979. Such a 
change would not affect the rain-gauge data, however. 

7. Discussion 

Daily rainfall over the SE U.S. (area B in Figure 3b) 
appears to have a significant weekly cycle; its amplitude 
r7 is about 15% of the mean. Statistical tests that take 
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account of natural variability with multi- day time scales 
indicate that the weekly cycles in both morning and after- 
noon rainfall and in atmospheric winds over the SE U.S. 
are extremely unlikely to be accidents of natural variabil- 
ity. Figure 4a shows that average rain rate over the SE 
U.S. is largest during the middle of the week and drops 
to its lowest value on Sunday — behavior very similar to 
anthropogenic pollution. 

Correlation does not prove causation, however. There 
are other types of human activity not directly contribut- 
ing to air pollution that vary on a weekly basis and 
that might affect storm dynamics. The heat produced 
in powering motor vehicles, industrial activity, and gen- 
erating electricity is too small to affect circulation on this 
scale. Jet contrails may be more numerous on weekdays, 
but would tend to reduce afternoon thermal instability 
rather than increase it. Irrigation of fields might have a 
weekly cycle, but it seems unlikely that this would change 
the available water in the atmosphere enough to account 
for the large changes in afternoon rainfall over the SE 
U.S. Over the drier SW U.S., where irrigation practices 
are more likely to affect atmospheric moisture content, 
a weekly cycle is not detectable. Other possibilities are 
discussed by de F. Forster and Solomon [2003]. 

These considerations suggest that summertime (mostly 
convective) rain amounts are increased, at least on large 
scales, by the increase in the kinds of pollution that vary 
with the day of the week because of human activity. The 
direct radiative effects of aerosols are unlikely to be re- 
sponsible, since aerosols tend to reduce surface solar heat- 
ing [Ramanathan et al , 2001; Koren et a/., 2004]; their ra- 
diative effects would tend to decrease midweek afternoon 
rainfall, in contrast with the observations. It is possible 
that concentrations of radiatively absorptive aerosols like 
black carbon could alter the vertical distribution of heat- 
ing and help initiate elevated convection in extreme cases 
of heavy black smoke [Rudich et al. , 2003]. Such clouds 
would be typically weaker than those starting from the 
boundary layer, however, and Koren et al. [2004] showed 
that the suppression effect on the triggering of convec- 
tion dominates the heating effect. It is therefore difficult 
to see how the heating due to aerosol absorption would 
explain the observed changes in the distribution of storm 
heights seen in Section 4. This would, however, explain 
the mid- 20th-century reversal of the weekly cycle, when 
the absorptive properties of the aerosols probably began 
to weaken while the CCN number concentration proba- 
bly continued to increase. On balance, based on what is 
known about the effect of aerosols on droplet formation 
in developing storms and on recent modeling studies, we 
believe the most likely explanation of the changes we have 
observed during the last few decades is the one proposed 
in the Introduction. 

It is important at this point to recall how the TMI is 
used to obtain rain-rate estimates. Over land the esti- 
mates are based mostly on microwave radiance measured 
by the 85-GHz channel. Since the 85-GHz signal is largely 
determined by the size and amount of ice aloft, which in- 
creases with the intensification of convection in summer 
storms, this method of estimating rainfall at the surface 
generally works quite well [Kummerow et al, 2001], How- 
ever, the size and amount of ice aloft can increase when 
more air pollution is added to the clouds while producing 
the same surface rain intensity [Rosenfeld and Ulbrich, 
2003; Khain et al, 2005]. It is therefore possible that the 
weekly cycle in the TRMM rain estimates may be due to 
changes in the ice aloft that are not necessarily accom- 
panied by such large changes in rainfall amounts at the 
surface. We note, however, that the gauge data appear 
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to show a weekly cycle smaller but still comparable in 
size to the cycle in the TMI data. The dynamical impli- 
cations of changes in ice aloft for cloud structure would 
be important, however, even if surface rainfall were to be 
unaffected. 

The weekly cycle in the Atlantic off the east coast of 
the U.S. (area C) is remarkable in that it is almost ex- 
actly opposite in phase to the cycle over the nearby land, 
and quite strong; see, in particular, Figure 9. Though it 
is not possible with the amount of data at hand to pro- 
vide a definitive explanation for this, we offer a tentative 
interpretation based on the following observations: 

1) As mentioned above, pollution aerosols tend to re- 
duce cloud droplet size, suppress early warm rainout and 
enhance the subsequent mixed-phase precipitation. Sea- 
salt aerosols restore much of the warm rain in polluted 
clouds and hence weaken the pollution effects over ocean, 
but do not reverse them [ Rosenfeld et al . , 2002]. 

2) Over water, the TMI senses rain much more directly 
than over land, and is less sensitive to ice aloft. Suppres- 
sion of the conversion of cloud drops to rain drops should 
not affect the TMI signal much because it is dominated 
by the amount of cloud and rain water combined and not 
very sensitive to their partitioning. 

These two considerations make it difficult to explain 
the observed midweek diminution of TMI rainfall over 
area C as a direct consequence of weekly variations in 
pollution in the area. Furthermore, we note that the 
weekly signal over land appears to die out near the coast 
and to reverse sign and become very strong over the open 
water. The phase does not seem to drift significantly to 
later times as one moves further eastward, which might 
have been expected if the suppression is due to pollution 
emanating from the continent. (CB report that they saw 
a steadily increasing delay in peak rainfall in the oceanic 
areas further and further east of area C, based on their 
all-season statistics. The almost simultaneous peak we 
see in the rainfall over area D in Figure 5a does not seem 
to be consistent with this, however. We examined the 
rainfall over an area east of D and adjacent to it: it has 
a slight weekly cycle, but its phase is not statistically 
distinguishable from that of the two oceanic areas C and 
D.) 

These facts suggest to us that what we are in fact 
seeing is dynamical suppression of midweek precipitation 
over the ocean in response to the invigorated convection 
over the land, as air pumped into the upper troposphere 
over the continent descends over the surrounding oceanic 
area, suppressing convection there and reducing its cov- 
erage and intensity. This is, in essence, a modulation 
of the monsoon-like land-ocean circulation in response 
to the aerosol-induced changes in convection. Validat- 
ing our proposed explanation will require quantitative 
model simulations of how changes in convection over the 
continent might affect offshore convection, and satellite 
documentation of the weekly cycle of pollution over the 
ocean. 

Baumer and Vogel [2007] point out that the weekly cy- 
cles they see in various variables do not seem to correlate 
well with likely local aerosol changes (they do not provide 
data for actual aerosol levels) inferred from the proxim- 
ity of urban areas to the stations they studied, and sug- 
gest that what they might be seeing is a regional rather 
than local effect. (It should be noted that they studied 
all-season averages, whereas we have presented statistics 
only for the summer season, and the physical mechanism 
proposed here is likely only to be applicable to summer 
convection. It is therefore difficult to make detailed com- 
parisons of their results and ours.) It is perhaps worth 
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noting that if the long-range effects of a weekly cycle 
in one locality can induce very different weekly cycles in 
nearby areas, as we have suggested may be occurring over 
the coastal Atlantic near the U.S., then local weekly cy- 
cles in precipitation may not be governed by local changes 
in pollution. 


8. Conclusions and Questions 

Based on summertime TRMM satellite, gauge, and 
reanalysis data, there appears to be strong statistical ev- 
idence that 

1) afternoon rainstorm activity is more vigorous dur- 
ing the middle of the week than on weekends over a sub- 
stantial portion of the southeast U.S., and the change is 
statistically highly significant; 

2) both the average area and intensity of rain increase 
in the middle of the week; 

3) there is a strong tendency for this weekly variation 
to show up in afternoon data in most (but not all) sum- 
mers; 

4) during the middle of the week, rain area is greater 
and hydrometeors tend to reach higher altitudes than on 
weekends, according to the TRMM radar; 

5) wind patterns change from midweek to weekend in 
a way that is physically consistent with the changes in 
convection implied by the satellite rain estimates; 

6) both gauge and reanalysis rainfall data are consis- 
tent with the satellite observations; 

7) the effect, which clearly must be anthropogenic, ex- 
tends over the nearby Atlantic, and is almost as strong 
there; and 

8) little evidence oTa weekly cycle is seen in the drier 
western half of the U.S. 

Based on the substantial amount of research docu- 
menting the influence of aerosols on cloud development 
and the weekly variations in aerosol concentrations, this 
evidence strongly suggests that air pollution invigorates 
storms in areas with large vertical instability, such as oc- 
curs over land in the summer, when there is an ample sup- 
ply of moisture. The enhanced convection is sufficiently 
vigorous that it may suppress convection in neighboring 
oceanic areas. The observations are consistent with the 
dynamical picture presented by Rosenfeld in Williams 
et al. [2002] and Andreae et al. [2004] as well as with re- 
cent model studies. The weekly cycle in rainfall thus 
seems to provide an opportunity for investigating the 
larger-scale climate implications of this picture. 

Many questions are" suggested by these results: 

1) Which constituents of air pollution are responsible 
for the rainstorm invigoration? 

2) Does the invigoration occur more around large ur- 
ban areas or over rural areas? [Wang, 2005, finds that 
the aerosol effect on model storm development tends to 
saturate at high concentrations, suggesting that storm 
intensification might be most noticeable in rural areas.] 

3) Is there a weekly cycle in rainfall in areas of the 
U.S. not visible to TRMM? Are there other areas in the 
world where a weekly cycle in rainfall would be expected? 

[Many areas do not seem to have strong weekly cycles of 
pollution. See Beirle et al . , 2003, for example. More- 
over, the absorbing aerosols in the most polluted regions 
of the world might neutralize or even reverse the CCN 
invigoration effect.] 

4) Is some part of the intensification of rain events 
during the past 30 years reported by Groisman et al. 

[2004] due to decadal changes in aerosol absorption and 
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CCN activity rather than the climatic effects of green- 
house gases? If weekly changes in human activity affect 
storm behavior, decadal trends in the same activities are 
likely to produce corresponding trends in storm behavior. 

5) Since air pollution is scavenged by rain, and rain 
appears to be increased by pollution (in the environments 
dealt with here) , the weekly cycle we observe is a complex 
product of these interactions and emissions. Can atmo- 
spheric models that incorporate aerosol effects reproduce 
the observed weekly cycle? 

6 ) Can atmospheric models reproduce (and help us 
understand) the day /night differences in the weekly cycle 
observed over the SE U.S.? 

7) More generally, how do aerosols affect the diurnal 
cycle of rainfall? 

8 ) Can atmospheric models reproduce the seesaw rela- 
tionship of the weekly cycle over land and over the nearby 
ocean? 

9) If there is a “weekend effect” on storm development, 
how much error is made in forecasts of precipitation and 
severe weather if forecast models do not include aerosol 
effects? 


Appendix A: Estimation of Statistical 
Significance of Weekly Cycles 

The likelihood that a weekly cycle seen in the data 
is an accident of natural variations in rainfall, which, of 
course, are very unlikely to favor any day of the week, 
has been estimated in two ways, both of which proved 
to give similar values for the statistical significance level 
of the observed cycles: The first method is based on one 
described in Bell and Reid [1993], and involves estimating 
the sampling error in a sinusoidal fit to the weekly cycle 
specified in Eq. (1). The time series is broken into 7-day 
chunks, each of which is fit to the linear version of Eq. 

( 1 ) with 3 unknown amplitudes, 

r(t) = ro + C 7 cos(u> 7 t) + sj sin(o; 7 t) (Al) 

with Tj = C 7 T Sj and <j ) 7 = (7/27r) aictan(s7 / C7) . If n 
weeks of data are available and provide n estimates of the 
coefficients C 7 and S 7 , then the error variance in C 7 and S7 
is estimated as the variance of the n estimates divided by 
n, assuming that the amplitudes are not very correlated 
from week to week. Under the null hypothesis 7*7 = 0, 77 
is distributed for large n as a chi-squared variable with 2 
degrees of freedom, and the probability that 7*7 exceeds 
R by accident is given by 

Prob(r 7 > R) = exp (-R 2 /a%) (A2) 

with erf = [var(c 7 ) + var(s 7 )]/n. As pointed out, for 
example, by Langmuir [1953] (p. 89) or Collier and Bow- 
man [2004] , if n is not large, it is more appropriate to use 
the Fisher F probability distribution, which takes into ac- 
count the fact that a 7 is estimated from a finite number 
of samples (n weeks), and gives, for this particular case, 

Prob(r 7 > R) = [1 + (R 2 /a 2 7 )/n']- n ' (A3) 

with n’ = n — 1. 

The second method we use is a bootstrap (resampling) 
approach that generates artificial (resampled) time series 
from randomly chosen chunks of data while attempting 
to preserve the time correlation in the data [e.g., Wilks , 

1997]. Resampling was carried out by using the original 
time sequence of observations, dividing it into chunks of 
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lengths 4, 5, or 6 days (randomly chosen for each chunk) 
and replacing each chunk with another of the same length 
randomly selected from those segments of the data whose 
local times of observation are within 2 hours of the chunk 
being replaced. This preserves the effects of the diurnal 
cycle of rainfall in the statistics, but scrambles any asso- 
ciations with the day of the week. A set of 10,000 syn- 
thetic time series and the associated fit parameters 
and <f> [ a \ a = 1, . . . , 10, 000, were obtained. The fraction 
of fits larger than the observed amplitude rj for the 
actual data was used to estimate the probability p that 
the observed amplitude could occur by accident, under 
the null hypothesis that there is no real weekly cycle. 
The value of <77 when resampling is used to estimate p is 
assigned the value 

g 7 = rj{— 1 np ) -1 / 2 . (A4) 

The variances of daily averages across these synthetic se- 
ries were also used to estimate the error bars for area B 
in Figures 4 and 7. 

The error bars in Figure 10c for the ratio of cumu- 
lative distributions of storm heights were estimated as- 
suming that the counts n(a) of how many storm heights 
exceed altitude a are counts of individual, independent 
events. If the total count of storm height values is, by 
definition, N = n( 0), then the cumulative distributions 
d(a) whose ratios are plotted in Figure 10c are given by 
Ci(a) = rii(a)/Ni, where i denotes either weekend (Sat- 
Mon) or midweek (Tue-Thu), and the ratios plotted are 
ctwt (a)/cssM (cl) . The assumption that the counts rii(a) 
are independent is certainly not true for low altitudes, but 
at high altitudes the counts are relatively small and the 
assumption is more reasonable. Denote expectations by 
angular brackets and write the observed counts m(a) as 

rn(a) = {rii(a)) 4- 5m(a) , (A5) 

where 5rii(a ) is the deviation from the expected mean 
due to the finite sample size. Because the counts are 
independent, n* (we omit specifying a henceforward) is 
Poisson distributed, with variance 


var(rii) = var(Jn*) = ( n* ). (A6) 

Then we can approximate sampling error in the ratio as 

£(ctwt/cssm) ~ (cssm) -1 £ctwt 

— (ctwt)(cssm) -2 <5cssm j (A7) 

from which it follows that the error variance is given 
by 

var(cTWT/cssM) ~ (cssm) 2 var(^CTWT) 4 - (ctwt) 2 

x(cssm) -4 var(^cssM). (A8) 

We have neglected the covariance term (£cssm<5ctwt) be- 
cause the cumulatives are from different days of the week. 

The variance of the cumulative-distribution errors Sci 
can be estimated using 

6ci « (Ni)~ 1 Sni — {rii){Ni}~ 2 5Ni (A9) 

and 

var(fe) = (( 5ci ) 2 ). (A10) 

Substitute Eq. A9 into A10. Use var(^A^) = ( Ni ) and 
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(SriiSNi) = (( Srii ) 2 ) = (n*), based on the Poisson as- 
sumption and the fact that Ni includes the counts in rii 
by definition. The variance of d is then given by 

var(<5cj) = cj {nj 1 - Nf 1 ) , (All) 

where we have replaced the expectations (rii) by their 
best estimates rii. Substituting this into Eq. A8 gives us 
our estimate of the standard deviation of the sampling 
error in the ratio, [var(cTWT/cssM)] 1/2 . 

Appendix B: TRMM Sampling Issues 

Because TRMM visits a grid box several days in a row 
within the same hour of the day, then visits the box for 
several days about one hour earlier, etc. [e.g., Negri et al . , 
2002], there is some danger that the satellite might, for 
example, have a tendency to overfly an area on weekends 
in the afternoon while visiting the area during the middle 
of the week in the morning. In other words, the diurnal 
variation of rainfall might induce a spurious weekly cycle 
in the TRMM statistics. This possibility has been ruled 
out by examining the number of TRMM observations of 
the areas in Figure 3b as a function of both the day of 
the week and the hour of the day: all hours of the day are 
about equally observed each day of the week. We have 
also tried replacing the actual time series of area-averaged 
rain rate with an artificial diurnal cycle characteristic of 
the area and found that the spurious weekly cycle in- 
duced by the diurnal cycle is two orders of magnitude 
smaller than the observed weekly cycle. Note, too, that 
an orbit- induced weekly cycle ought to have shown up in 
the southwest- U.S. area- A statistics, and it does not. 

TRMM’s orbital altitude was increased from 350 km 
to 402.5 km in August 2001 to conserve its remaining fuel. 
Changes in rain statistics after the boost are subtle, and 
unlikely to affect the amplitude or phase of the weekly 
cycles seen in the data. 

TRMM’s orbital plane is inclined 35° with respect to 
the Equator, and as a consequence the TMI and PR in- 
struments view, respectively, latitudes around 31° and 
34° roughly 3 times as often as they do Equatorial lati- 
tudes [e.g., Bell et al . , 1990]. TMI coverage extends only 
up to about 38.4° latitude, and the PR only up to about 
36° (prior to TRMM’s orbital boost). Spatial averages 
of TMI and PR data used here are numerical averages 
over individual fields of view of the instruments, and so 
are weighted most heavily towards the best observed lati- 
tudes, unlike the spatial averages of gauge and reanalysis 
data. This choice in averaging was guided by the fact that 
if a weekly cycle in rainfall was present at equal levels ev- 
erywhere in an area, then frequency-weighted averaging 
would produce the highest signal-to-noise ratio. This dif- 
ference in weighting of the satellite and the gauge and 
reanalysis data might in principle make comparisons of 
the averages difficult, but in an experiment we conducted 
in which the rain-gauge data were averaged with weights 
proportional to the TMI sampling frequency we found 
little change in the results when compared to those with 
equal weighting of all gauges. 
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Phase, Significance Level of Weekly Cycle of PM 10 (JJA, 1998-05) 



Significance Level p 

Figure 1 . Day of week of maximum of sinusoidal fit to daily PMio measurements at 132 sites in the 
U.S., and to daily averages of hourly PMio measurements at 254 sites. Both color and lengths of arrows 
are used to indicate the statistical significance level p of the weekly cycle. Lengths are proportional to 
(— Inp) 1 / 2 . Keys to directions are shown in clock format at right, and to length scales by arrows to right 
of color bar. 
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Phase, Significance Level of Weekly Cycle of PM 2 5 (JJA, 1998-05) 



Significance Level p 

Figure 2. Day of week of maximum of sinusoidal fit to daily PM 2.5 measurements at 212 sites in the 
U.S. Both color and lengths of arrows are used to indicate the statistical significance level p of the weekly 
cycle. See Figure 1 caption for explanations. 
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a) Phase of Weekly Cycle (J JA, 1 998-2005, v6) 



Sat Sun Mon T ue Wed Thu Fri 


b) Averaging Regions 



Figure 3. Study region, a) Gridded phases of weekly cycle. Colors for each grid box indicate day of 
the week with maximum rainfall, based on fits to a 7-d sinusoid of averages of TRMM rain estimates for 
8 summers. Grayness of colors is based on the significance level p for the amplitude, with colors along 
top of color bar indicating significance level p < 0.05. Areas with rain rates less than 0.01 mmh~ 1 axe 
masked: b) Five averaging areas chosen to improve the signal-to-noise ratio for cycle detection. 
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Sat Sun Mon Tue Wed Thu Fri Sat Sun Mon Tue Wed Thu Fri 



Figure 4. Rain statistics (JJA, 1998-2005) from TRMM 
for each day of the week for areas A-C delineated in Fig- 
ure 3b. Area-B error bars are 1-sigma confidence lim- 
its determined using bootstrap method, a) Average rain 
rate, b) Fraction of TMI footprints with rain, c) Ratio 
of results in panels a and b, indicating intensity of rain 
where it rains (i.e., conditional on R > 0); note vertical 
scale in panel c. 
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Figure 5. Rain statistics as in Figure 4 for axeas C-E in 
Figure 3b. Statistics for area C in Figure 4 repeated here 
to aid comparisons. Note vertical scale change in panel 
c. 
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Weekly Cycle Phase and Significance Level p 
Sun 



C - Atl, Coastal 
D - Atl, East of C 
E - Gulf of Mexico 


Figure 6. Phase and significance level of weekly cycles of 
mean rain rate for 5 areas in Figure 3b. Angular position 
is determined by day of week with maximum rain rate, 
based on sinusoidal fits to plots in Figures 4a and 5a. 
Radial distance is determined by significance level p of 
weekly cycle and is proportional to (— lnp) 1/2 . 



Figure 7. Half- day- average SE-U.S. (area-B) rain rate 
for mornings (0000-1200 LT) and afternoons (1200-2400 
LT). Error bars are 1-sigma confidence limits estimated 
using resampling. 
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Annual Phase, Significance of 
a) TMI Weekly Cycle 

Sun 



JJA, 1998-2006, Afternoon 


b) TRMM 3B42 Rain Rate 

Sun 



JJA, 1 998-2006, Afternoon 


Figure 8. Phases and significance levels of weekly-cycle 
fits to average afternoon (1200-2400 LT) rain rate over 
SE U.S. (area B in Figure 3b) for 9 individual summers 
for 1998-2006. a) TMI v. 6 data alone, b) 3-hourly 
TRMM merged satellite product (TMPA). 
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Figure 9. Average rain rate as a function of the day of 
the week and hour of the day (local time) using data from 
Jun-Aug, 1998-2005.. Hourly values have been smoothed 
with a “centered” running 4-h average, r* = [(l/4)(rt-2+ 
n-i+n+n+i)+(l/4)(rt_i+rt+rt+i+n+2)]/2. Because 
the day advances by 1 with a 24-h shift, data for hours 
00-23 have been repeated for horns 24-47 but assigned 
to the next day of the week. The right-hand vertical axis 
is therefore shifted by 1 day relative to the left-hand axis, 
a) Averages for area B in Figure 3b. b) Averages for area 
C. 
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Morning Storm Height Distributions 



Frequency / km 

Afternoon Storm Height Distributions 



Frequency / km 

Ratio (Midweek / Weekend) of Probabilities 



Figure 10. Frequency statistics of PR storm heights 
over area B using 1-km- altitude bins, derived from 
TRMM PR product 2A23, v. 6. a) Morning (0000-1200 
LT) storm heights. Circular dots show the distribution 
of “midweek” storm heights, while triangles show “week- 
end” distribution. Dashed curve shows the difference, 
b) Afternoon (1200-2400 LT) storm heights, c) Ratio 
of midweek cumulative distributions to weekend cumu- 
lative distributions, representing the relative probability 
that midweek storm heights exceed weekend heights, with 
two-sigma error bars. (Cumulative distributions defined 
as integrals of the probability distributions from the high- 
est altitude downward.) 
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TMI Precip (mm / day) — • - 300 hPa Div (1 O' 6 s -1 ) 
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Figure 11. Wind-field statistics (summer, 1998-2005) 
from NCEP2 reanalysis over areas B and C in Figure 3b 
as a function of the day of the week. Average wind con- 
vergence over the area at 1000 hPa indicates net inflow 
into the area near the surface. Average divergence at 300 
hPa 9 km) indicates net export of air out of the area 
at altitudes typical of storm tops. Mid- atmosphere (~ 5 
km) average vertical velocity at 500 hPa is represented 
by — u>j the time rate of change of air pressure of an air 
parcel (sign reversed). Daily TMI anomalies, based on 
Figure 4a, are superimposed for comparison. Differences 
from time means are plotted for each quantity. Statisti- 
cal significance levels for the 3 wind fields for area B are 
p = 0.024, 0.045, and 0.011 y respectively; for area C they 
arep = 0.047, 0.41, and 0.20. 
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-TRMM (TMI) Precip 
-Rain Gauge (GHCN) 


■NCEP-1 Precip 
-NCEP-2 Precip 
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Figure 12. Area-averaged rain rates from TMI, NCEP 
reanalyses versions 1 and R-2, and GHCN rain gauges for 
each day of the week. 1-sigma error bars are estimated 
using resampling for TMI data and (7 / 4) 1 1 2 <t 7 for gauge 
and NCEP data, where cr 7 is obtained as described in 
Appendix A. a) Area-B averages, b) Area-C averages. 
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Figure 13. Statistics of weekly-cycle fits to area-B- 
averaged GHCN rain-gauge data. Signal-to-noise ratios 
7*7 /cr 7 are plotted in the top panel (see text for explana- 
tion of horizontal lines labeled with p), and phases (j ) 7 in 
the bottom panel, a) Ratio of t^/c^ from fits to 8 con- 
secutive summers of data ending on the year plotted, b) 
Corresponding phase (j ) 7 of fits. 
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Figure 14. As in Figure 13, but showing statistics of weekly-cycle fits to area-B- averaged NCEP1- 
reanalysis precipitation (panels a,b), 1000- hPa horizontal- wind convergence (panels c,d), 300-hPa 
horizontal- wind divergence (panels e,f), and 500-hPa vertical- wind (— a;) (panels g,h). 
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Table 1. Parameters of weekly-cycle fits to area-averaged rainfall 


Region 

r 0 a 

r 7 a 

h b 

P 

Area A (SW U.S.) 

0.129 

0.0035 

5.4 

0.86 

Area B (SE U.S.) 

0.174 

0.025 

4.7 

0.10 

Area C (Coastal Atl.) 

0.221 

0.042 

0.9 

0.035 

Area D (Atl. East of C) 

0.165 

0.028 

0.6 

0.10 

Area E (Gulf of Mex.) 

0.098 

0.016 

1.1 

0.18 


a Units: mm h _1 . 

b Time of the week when sinusoidal fit peaks, with 0 < $7 < 7 and 0.0 — > 0000 Saturday. 
Table 2. Average fraction of area with rain seen by PR 



Morning 

Afternoon 

Sat-Mon 

Tue-Thu 

0.0297 

0.0245 

0.0415 

0.0579 


Fraction of PR footprints containing rain with high confidence, from TRMM PR product 2A23, for morning (0000-1200 LT) and 
afternoon (1200-2400 LT) hours for weekend and midweek periods, over area B for JJA 1998-2005. 



Final Year of 1 5-Year Window 

Figure 15. Signal-to-noise ratios r 7 / 0-7 of weekly- cycle fits to area-B- averaged summer GHCN gauge 
data and NCEP1 reanalysis estimates of precipitation, surface wind convergence, 500-hPa vertical wind 
velocity, and 300-hPa wind divergence using a 15-year window. Value plotted for year y represents fit to 
years y — 14 to y. 
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Table 3. Statistics of fits to rain data (1998-2005) 



TMI 

Precip 

Rain 

Gauges 

NCEP1 

Precip 

NCEP2 

Precip 

ro 

4.18 

Area B 
3.31 5.74 

5.21 

r 7 

0.60 

0.35 

0.46 

0.47 

cr 7 

0.40 

0.22 

0.26 

0.32 

(j> 7 

4.7 

5.6 

4.9 

5.5 

P 

0.101 

0.089 

0.044 

0.109 

ro 

5.31 

Area 

C 

4.16 

6.04 

r 7 

1.06 

- 

0.37 

0.93 

<J 7 

0.55 

- 

0.26 

0.44 

<j> 7 

0.9 


0.5 

0.9 

P 

0.035 

- 

0.124 

0.011 


Parameters of fits to JJA daily rain-rate data from TMI, GHCN gauges, NCEP1 and NCEP2 reanalyses. Parameters r 0 , r 7 , and 
(77 have units mm/day. The sinusoidal fit peaks on day (j> 7 , with 0 < <j> 7 < 7 and 0.0 —► 0000 Saturday. 



